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O ■ Abstract 



Non-equilibrium Coulomb effects in resonant tunnelling through deep impurity 
states are analyzed. It is shown that Coulomb vertex corrections to the tunnelling 
transfer amplitude lead to power law singularity in current- voltage characteristics. 
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Localized states of individual impurity atoms and interacting impurity clusters can 
play the key role in tunnelling processes in small size junctions and often determine the 
behavior of tunnelling characteristics in STM/STS contacts. Now it is evident that in tun- 
\ nelling junctions of nanometer scale there exists non-equilibrium distribution of tunnelling 
electrons which changes local density of states and tunnelling conductivity spectra. Some 
^ ■ interesting effects, such as resonance structure of tunnelling conductivity inside semicon- 
ductor band gap, increased value of observed band gap and non-equilibrium interaction 
of neighboring impurity atoms have been recently investigated experimentally and the- 
oretically analyzed [p], ||, H, U- But all these effects are caused by local changes of the 
initial density of states in the contact area due to interactions of non-equilibrium parti- 
cles. The modification of tunnelling amplitude by the Coulomb interaction of conduction 
electrons in metallic tip with non-equilibrium localized charges was ignored. It is shown 
in the present paper that corrections to the tunnelling vertex caused by the Coulomb 
potential can also result in nontrivial behavior of tunnelling characteristics and should 
be taken into account. One encounters with effects similar to the Mahan edge singu- 
larities in the problem of X - ray absorption spectra in metals JS|. The effect is well 
pronounced if tunnelling rate from a deep impurity level to metallic tip % is much larger 
than relaxation rate 7 of non-equilibrium electron distribution at localized state. This 
condition can be realized experimentally for a deep impurity state in the semiconductor 
gap. Direct tunnelling from such states to semiconductor continuum states is strongly 
reduced due to wide barrier formed by surface band bending. Relaxation rate connected 
with electron-phonon interaction can be estimated to be of the order 10 8 — 10 10 l/s at 
low temperatures As to 7^ — it is parameter, which can be varied in STM/STS 
experiments changing tip-sample separation. Since tip-sample separation is comparable 
with atomic scale, jt often exceeds the relaxation rate for deep impurity states. Typical 
experimental value of tunnelling current 1 nA corresponds to jt ~ 10 11 — 10 12 l/s ||. 
As it will be shown below for 7t 3> 7 the impurity level becomes nearly empty when the 
value of applied bias voltage approaches the impurity energy. So the core hole Coulomb 
potential is suddenly switched on and the tunnelling amplitude is changed. One might 
expect in this situation a power-law singularity in current- volt age characteristics near the 
threshold voltage. 
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The system semiconductor - impurity state - metallic tip can be described by the 
Hamiltonian H: 

H = H L + H R + H imp + H T + H int (1) 

where: 

H R = Y ( £ k ~ V) C k*Cka, H L = J2( £ p-^~ eV )4o°P° ( 2 ) 
ka po- 

describes the electron states in the metallic tip and the semiconductor correspondingly, 
c tA c ka) and Cp a (c pa ) describe creation (annihilation) of electron in states (ka) and (per) 
in each bank of the contact. 

H imp = £ d C da C da + U n d(j n d _ a (3) 
do 

corresponds to a localized impurity state. We consider "one electron neutral impurity" - 
the impurity level is single occupied at zero applied voltage due to the on-site Coulomb 
interaction. But in the case of large tunnelling rate to the metallic tip the on-site Coulomb 
repulsion of localized electrons can be omitted, if we analyze the behavior of tunnelling 
current at applied voltage close to the impurity energy s d . Because in this situation the 
impurity state becomes nearly empty above the threshold value of the applied bias. Let 
us also point out, that the Kondo regime is destroyed on Anderson impurity for the values 
of applied bias near the threshold |7|, ||. In this case the Kondo-effect is not responsible 
for any unusual features of the tunnelling characteristics . 

Tunnelling transitions from the impurity state to the semiconductor and the metal are 
described by the part: 

H T = (T kd c\ a c du + T pd c+ a c da ) + h.c. (4) 

kp 

And, finally, part H int includes the Coulomb interaction of the core (impurity) hole with 
conduction electrons in the metal. 

H int = Y w kk>ct a c k < a (l - c da ,c da/ ) (5) 

kk'aa' 

Hamiltonian Hi nt appears as many particle interaction and describes rearrangement of 
conduction electrons in the potential of the hole, suddenly switched on by tunnelling 
transition of the impurity electron. Since we are far from the Kondo regime it is suffi- 
cient to consider tunnelling current in the lowest order in the tunnelling amplitude T^d- 
Scattering by the impurity hole Coulomb potential does not change electron spin. Thus 
in the lowest order in T kd we can consider renormalization of the tunnelling amplitude 
independently for each spin - the same one for conduction and impurity electrons. It is 
also reasonable to use for simlicity an averaged value of screened Coulomb interaction 
describing s - wave scattering of conduction electrons by a deep hole Wkk' = W. 

Edge singularities in the tunnelling current can be analyzed by means of diagram tech- 
nique for non-equilibrium processes. Using Keldysh functions G < the tunnelling current 
can be determined as (we set charge e = 1): 

I(V) = Im(J(V)) , J(V) duT kd Gl< (6) 

k,a 
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where we have defined tunnelling "response function" J(V). If the Coulomb interaction 
is neglected one can obtain the usual expression for this response function in the lowest 
order in T kd : 

J°(V) =«£/ duTUGl<G^ + GffGS?) (7) 



k,cr 



Substituting the correspondent expressions for the Keldysh functions M and performing 
integration over k we get: 



J°(V) = lt / du 



n° k (u) n d (u)(-i(j + j t )) 



oo + eV -e d + i(7 + 7*) (u + eV - £ d ) 2 + (7 + 7t) s 

where the tunnelling rate jt — T kd y, and z/ is unperturbed density of states in the metal- 
lic tip. Kinetic parameter 7 corresponds to relaxation rate of electron distribution at 
the localized state. In the suggested microscopic picture (eq.(4)) this relaxation rate 
is determined by small enough electron tunnelling transitions from the impurity to the 
semiconductor continuum states 7 = T pd v p . (In general 7 can include different types of 
relaxation processes.) 

Non-equilibrium impurity filling numbers n d (u) are determined from kinetic equations 
for the Keldysh functions G < : 

jn°{u)+j t n° k {uj) 

n d {uj) = 9 

7 + 7t 

As it was explained in the introduction for a deep impurity level the relation 7$ >> 7 is 
quite possible. Then n° p {e d ) = 1, while n d (e d ) << 1 and there is really a core hole in the 
impurity state. Thus for low temperatures one can obtain from Eq.@: 

J°(V) = lt ln{\X\)+i^- [arcctg((eV - e d )/( 7 + lt )) - arcctg((-e d ) / '(7 + 7*))] (10) 

7 + It 

where 

X = {eV-e d + i{j + lt ))/D) (11) 

and D is the band width for electrons in metal. 

The usual form of the tunnelling current is of course reproduced from Eqs.(|6|.|8| JT0|) . 

Now let us consider renormalization of the tunnelling amplitude and vertex corrections 
to the tunnelling current caused by the Coulomb interaction between the impurity core 
hole and electrons in the metal. Many particle picture strongly differs from the single- 
particle one near the threshold voltage. First order corrections due to the Coulomb 
interaction (the first graph in Fig. 1(a)) has logarithmic divergency at the threshold 
voltage eV = e d , which is cut off by the finite relaxation and tunnelling rates. 

J\V)=tJ2 J dujT kd ( - G u kk G dd T kd + + Gl*G°<T l kd -) (12) 

Tunnelling matrix elements are changed by the Coulomb interaction: 

Tlr = E / du,T kd W(G%<G°J + Gl*G°<) (13) 

fc,(T 
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Figure 1: Coulomb corrections to T kd . Solid lines represent G k and dashed lines — G d - 
a) Ladder approximation, b) Parquet graphs (Coulomb wavy lines are overdrawn as black 
circle vertexes) 



If we look at Eqs. (|7],§), h becomes clear, that logarithmic contribution comes from 
the first combination of the Green functions: G kk G dd . In what follows we retain only 
logarithmically large parts, assuming that \ln(( , ~f + j t ))/D) \ ^> 1, so only these combina- 
tions of the Green functions are the most important in perturbation series. Then from 
(13) we obtain that tunnelling amplitude contains logarithmic correction: T kd = — T kd L, 



Tkd + = ~T kd , where factor L: 



L = (Wu)ln(X) (14) 



In high orders of perturbation expansion ladder graphs (Fig.l a) are the simplest 
"maximally singular" graphs. But this is not the only relevant kind of graphs. If we 
look at the first graph in Fig. lb, we notice, that a new type of "bubble" appears, which 
is logarithmically large for small "total" energy (u + Ui). The important point is that 
relevant region of integration over u and uj\ is region of small uj. It is just this region 
which gives essential contribution to logarithmic factor L in any other pair of G kk G dd R \ 

It means that the central bubble also contributes an additional logarithmic factor to 
the total result. In this situation, which is not new in physics, one should retain in the n - 
th order of perturbation expansion the most divergent terms proportional to {Wv) n L n+1 
. For the first time such method was developed by Dyatlov et al |J. It was shown 
that for a proper treatment of this problem one should write down integral equations for 
the so-called parquet graphs (Fig. lb), which are constructed by successive substitution 
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the simple Coulomb vertex for the two types of bubbles in perturbation series. These 
equations represent some extension of the ordinary Bethe-Salpeter equation and describe 
multiple scattering of conduction electrons by the core hole Coulomb potential in the 
two "most singular" channels. The integral equations can be solved with logarithmic 
accuracy, as it was done, for example, by Nozieres ]TD, |TTJ for edge singularities in X - 
ray absorption spectra in metals. 

The solution of the "parquet" equations contains nothing new in our problem. So we 
present the result without getting into technical details. Summing up the most divergent 
graphs with logarithmic accuracy one can obtain the following singular part of the response 



function 10 



Then the tunnelling current near the threshold voltage can be expressed as: 

It 



I(V) 



2Wv 



(eV -Edf + ilt + lf 



sin (2Wv<j>) (16) 



where = o,rcctg( e ^ + ^ d ). If we consider a deep impurity state in the gap of the 
semiconductor (below the Fermi level) and positive tip bias voltage, then: < 0, eV < 0. 
So, the phase is a step-like function varying approximatly from to tt, when the applied 
bias crosses the threshold eV = Ed- Since we retain only the most logarithmically large 
terms in the tunnelling current, Eq.(^) is valid only if \eV — Ed\ ^ D. In the absence 
of the Coulomb interaction (W = 0) this singular part reduces to the usual first order 
contribution arising from the first term in Eqs.(^.^|). 

In summary, if the tunnelling rate to the metallic tip exceeds the relaxation rate of 
localized electrons, the Coulomb interaction of the core hole and the conduction elec- 
trons in the metal strongly modifies the tunnelling transition amplitude and leads to: 
i)non-monotonic behavior of current voltage characteristics; ii)power-law singularity of 
the tunnelling current and conductivity when the value of the applied voltage approaches 
the impurity level energy; iii) current voltage characteristics can be rather asymmetric, 
because of different dependence of the phase factor on the applied bias below and above 
the threshold value. Power law singular behavior of the tunnelling current is sensitive to 
the values of the tunnelling and relaxation rates, as well as to the value of the Coulomb 
interaction W. So, different exponents in power dependencies of the tunnelling current 
on the applied voltage can appear with changing tip- sample separation. Some current- 
voltage characteristics obtained for typical values of parameters are shown in Fig. 2. 

It seems also possible to set up an experiment with negative impurity charge and 
negative tip voltage close to the value Ed + U. In this case W > 0, and the Coulomb 
corrections to the tunnelling amplitude result in power-low behavior of the tunnelling 
current with the opposite sign exponent in Eq. (p~6|) . The tunnelling current is suppressed 
near the threshold, compared to the noninteracting case. This behavior is shown in Fig. 
2b. 

Experimental STM/STS investigations of deep impurity levels on semiconductor sur- 
faces give evidence of the existence of the described effects. Some STM images demon- 
strate non-monotonic dependence of tunnelling current on applied bias voltage M . With 
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Figure 2: Current- voltage curves for typical values of dimensionless Coulomb and kinetic 
parameters. Current is measured in dimensionless units I/ej. a) w = Wv < 0,Ed = 
OAeV, b) w — Wv > 0,ed + U = OAeV, •jt/l = 3, Sd/l = 40. Dashed lines correspond to 
W = 0. Experimental STM image of Cr impurity on InAs (110) surface is shown in the 
inset for V=0, V=0.5(V), V=1.5(V) in sequence. 
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increasing of bias voltage impurity atom is "switched on" in STM image - appears as a 
bright spot. But further increase of the tunnelling bias "extinguishes" the brightness of 
the impurity atom and it is seen as a dark spot in the STM image (see inset in Fig.2). 
According to the present model it can be explained by decrease of the tunnelling current 
caused by the Coulomb vertex corrections to the tunnelling amplitude for some types of 
impurities. 

We thank L.V.Keldysh for helpful discussions. This work was supported by the RBRF 
grants N 00-15-96558, 00-02-17759 and program " Nanostructure" . 
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